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TECHNICAL NOTE 2471 1 

UNSTEADY LAMINAE BODNDAEY-IAYEE FLOW 
By Itaiiklin E. Moore 


SDMMAEY 

The leuainar compressihle hourudary layer oyer an inaulated flat 
plate moying -with a tlme-depencLent yeloolty has heen analyzed in detail. 
A group of parameters arise idiich, if large, proyide that the classical 
"starting from, rest" solution applies, and, if small, that the metion 
is q.uaBi-steady, These parameters relate to the time required for 
temporal changes to diffuse throu^ the boundary layer. Deyiations 
from the q.uasl-steady yeloolty and temperature profiles haye been 
computed. 

Unsteady laminar flows with pressure gradient and, probably, 
unsteady turbulent boundary layers are goyemed by similar parameters, 
which may be estimated In order to proyide a criterion as to whether 
q.uasi -steadiness may be assumed for a giyen problem. The case of a 
fluctuating yelocity field passing oyer a flat plate is discussed from 
this point of yiew. 


INTEOnjCTION 

The unsteady laminar boundary- layer flow oyer yarlous bodies 
starting from rest has been analyzed by seyeral inyestigators , the 
first of whom was Blaslus (see paragraph 65, reference l) . These 
treatments haye been concerned with such matters as the onset of 
separation on an airfoil and the transient effects of impulslye start; 
they employ an expansion of some suitable quantity in powers of the 
time elapsed since the start of motion. Thus, only the earliest phase 
of motion is considered. 

Unsteady boundary- layer flow for longer times elapsed since 
starting and for hi^er speeds req.uire consideration. Eor example, 
the flight speed of the usual rocket missile yaries continuously oyer 
the entire trajectory, and thus the Important boundary- layer effects 
of skin friction and heat transfer must, in principle, be regarded as 
unsteady for the entire fli^t. Eecently, Eaye (reference 2) has 
discussed the problem of unsteady heat-transfer effects associated 
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with nLissile flight, aaad has tacitly assumed that at hl^ speeds the 
■boundary layer responds with no time lag to changes in stream 
Telocity; that is, that the boundary layer at any Instant is that which 
would be associated with steady motion at the stream conditions pre- 
vailing at the same instant. Such a boundary layer is hereinafter 
called q.uasi-steady.“ - 

Another type of problem of-unsteady boundary- layer flow for high 
speed and long elapsed time conceims the response of a boundary layer 
to time -dependent velocity or pressure fluctuations in an otherwise 
steady outer potential flow. 

In the present report of research conducted at the HACA. Lewis 
laboratory, the case of compressible laminar flow over a semi-infinite 
flat plate" in rectilinear accelerated fll^t through still-air is 
treated in detail; the flight speed 1s considered to vary with time in 
a continuous but otherwise arbitrary manner. The solution to this 
problem, which is intended as an idealization of missile fll^t, will 
be interpreted to provide insist into a more general class of 
unsteady boundary-layer flows, both laminar and turbulent. 

The notation used in this report is described in appendix A. 


ARALIBIS OS’ ILO¥ O'VER FLAT PLATE IN ITNSTEADY jBLIGBT 
Eq.uatians of Motion 

The motion to be considered is that of a semi-infinite flat plate 
moving in a straight path, normal to its leading edge and in its own 
plane. Aside from the disturbance due to the motion of the plate, the 
medium (air) is at rest. The speed of flight may vary with time. A 
rectangular coordinate system (fig. 1(a)) is chosen to be at rest in 
the fluid, such that" the plate moves with a velocity U(t) in the 
negative x-direction, and y is measured normal to the plate. The 
origin of coordinates is taken to be the leading-edge location at 
t = 0. 

The Prandtl boundary-layer assumptions are presumed to be valid 
for this flow and, in particular, the pressure is taken to be constant 
throu^out the fluid. The Prandtl number and specific heats are 
assumed to be constant. Thus, the appropriate eq.uations of 
motion of a compressible fluid are; 
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P0 s= constant 



= 0 



These eq.uations may he written in a new coordinate system 
(fig. 1(h)) fixed with referance to the plate, with the origin at 
the leading edge. Temporarily denoting hy primes the new coordinates 
and the corresponding velocity components, the appropriate trans- 
formations are; 


V* 5 v; u* s u + Uj 

rt 


X’ E Z + 

Eq.uations (l) hecome; 




U dt; y‘ = y; t’ s t 



„ ro j. 

a L 

5u*\ 

6j' 


dy* V 


he 

1 s 


, _u_ 


“ Pr dy’ 


•4“ 






^ + ^ = ° 


J 


( 2 ) 


p6 = constant 


Hereinafter, the primes will not he used, and x, y, t, u, and v 
are understood to he measured relative to the coordinate system fixed 
in the plate (fig. 1(h)). 

The temperature-viscosity relationship discussed hy Chapman and 
Euhesin in reference 3 is employed; that is, it is assumed that the 
variation of viscosity with temperatur'e may he approximated hy the 
relation 
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(3) 


where C le‘ a constant obtained by matching equation (3) to the 
Sutherland fomula at (for example) the wall. This matching yields, 
according to reference 3, 


C = 


f^f Qq. + 216 ° S \ 

V ©00 +-216“ sy 


(4) 


The equations of motion are sinq)llfied by the use of a transformation 
of coordinates, similar to that employed in reference 4, and of a 
stream function, as followe: 


,7 


Y 3 


(p/Poo) cLy; X s x; T s t 




*^0 

P«x) Sijr ^ 

^ ^ T^ = '5y 



> (5) 



"Sc ^ 





Relations (3) and (S) transform the equations (2) of motion to yield: 


"^X^YI U>(T) t ViYYY 

©51 + ~ ^X®Y “ ^^®YY ^ (^YY^ 

The appropriate boundary conditions on if are 

\|fy(X,«>) = ifjiO,Y) = U(T) (7a) 

l|fy(X,0) = 0 (7b) 

y(X,0)=0 (7c) 


The boundary conditions on ©(X,Y) must be chosen with regard to the 
particular problem under consideration. (Subscripts and primes denote 
partial and ordinary differentiation, respectlTely. ) 
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Solutions for V, in the Case of Cooistant Acceleration 

If the Telocity is glTen hy the relation TJ(T) = AT, time T 
"being measured from rest, the eq.uatlons (6) may he -written in similarity 
form, defining 


t 3 fJT At 3/2 g(z ,0 

where 


Z = Y/sJW; {hX/AT^ 

These definitions enable eq.uations (6a) and (7) to he -written as 
follows : 


" I ^ZZ ~ ^^®Zt ^^z{ ~ ®£®zz “ ^ ®zzz 

g 2 (z,o) = = 1 ; 82(0^0 » g(o,C) - 0 ( 9 ) 

Eq.uation (8) is obtained hy considering u = — 0 to he a constant. 

Subject to the limitations imposed hy this assumption, which is dis- 
cussed in appendix B, this analysis may he considered applicable to 
compreselhle flow. 

Case t ^ ^ 1 » - If i * s/AT^» 1, the function g(Z, {) may he 
expanded in reciprocal powers of C as foUo-ws: 

g = go(Z) + Y gi(Z) + i ggCZ) + . . . 
go'(*) = gi'(*) = gg’(*) = . . . = 0 

go'(O) = 0; gi’(O) = g 2 *( 0 ) » . . . = 0 

ggCO) = Oj g^(0) = gg(0) = . . . = 0 

"When, in eq.uation (8), successlTe reciprocal powers of { are eq.uated, 
this definition and set of boundary conditions yield a sequence of 
differential equations, the first of which is 

so' ”i + so'" 
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and tlie solutions 


g, 


0 


. = _ i Z2 ^ + (Z2 + 2) 


pZ/2 


-z2 

e dz 


gn = go = Et; = • • • = 0 


This approach proyldes a solution equivalent to tha,t for the initial 
motion of a flat plate with constant acceleration presented in 
paragraph 66 ^ reference 1. 

This analysis clearly pertains to the very earliest stages of 
motion since the function g(Z, £) is expanded, in effect, in powers 
of time elapsed since the stasrt of motion. In this connection, it may 
he noted that the inertia teams in tjie momentum equation drop out as a 
consequence of assuming {»1. 

The boundary condition g2;(Z,0) = 1 has not been applied. Thus, 

no consideration is given to the growth of the boundary layer in the 
X-directlon, and the effect of the leading edge is ignored. The 
quantity { is equivalent to the ratio of the distance (X) aft of 
the leading edge to the distance, J i AT^j that the plate has 

traveled. (See fig. 1(a).) If {>>1, it is therefore proper to 
consider that the effect of the leading edge has not yet been felt at 
the station X. 

Case ^<<1 . - If, however, the distance traveled is much larger 
than distainoe aft of the leading edge, that is, if C<< 1, then 
account must be taken of the growth with X of the boundary layer. 

The appropriate eapansion of g(Z, {) for small values of { is 
therefore 



subject to the boundary conditions 

g0'(®“)=2; gi»(») 

go’(o)=o^ 

go(o) = Oj g^(0) = 

Substitution of this expansion of 
sequence of differential equations 


— S2 * ~ “ 0 

= gg*(0) = . • . = 0 

g2(0) — y • • ■ 0 

g(Z,C) into equation (S) yields a 
the- first of which is 


2229 



2229 


NACA W 2471 


7 


goSo“ + so’“ = 0 

!TMe differential eq.uatlon and the corresponding set of "boundary condi- 
tions are precisely those -which apply in the case of steady flo-w oyer 
a semi -infinite flat plate. (See reference 5.) In the steady case, 

the argument of gQ is 
argument is 

1 _Z_ 1 Y / / "Y" _ 1 _Y- /sn 

Therefore, at any instant when {<<1, the boundary- layer flo-w for 
constant acceleration is nearly that which would obtain in steady flow 
at a velocity eq.ua! to the Instantaneous velocity of the accelerated 
plate; as 0, the eq.ul-ralence becomes precise. The correction terms 
arising when C is small but not zero -will be considered in a subse- 
q.uent section. 


Y 


— In the present analysis, the 


Parameters for Unsteady Flat-Plate Flow 

Before the flow over a flat plate -with an arbitrary velocity-time 
relation is considered, it is desirable to determine the governing 
parameters. Fcpm the important physical q.uantltles Z, V, U, U’(T), 
U''(T), . . . U'°-^(T), . . . the following dimensionless q.uantlties 
may be constructed (in addition to Reynolds number and Mach number) : 

2U‘ . Z%*« Z°u(“-) 

^2 " -gS ' • • • " U ^+1 ^ 

Quantities (lO) and Reynolds number are therefore the q.uantitles eapjected 
to govern the unsteady Incompressible boundary layer on a flat plate with 
arbitrary (but differentiable) unsteady flight speed. These q.uantltles, 
in the same form, also govern the corresponding compressible flow, as 
will be shown in the next section. 

If U = AT, fes in the pre-vlous section, then the q.uantltles (lO) 
specialize to the single parameter This q.uantlty ^ has 

At2 

been noted to be eq.ui-mlent to the ratio of the distance aft of the 
leading edge to the distance traveled by the pla-be, 

A more direct physical meaning may be attached to the q.uan- 
titles (lO), and hence to the special case &b follows; 
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If some physical quaatity were suddenly to change at the "boundary 
of a fluid region of thidmess 8 in which conditions were initially 
uniform, this change would diffuse, with •time, throu^out the reglonj 
uniformity would ultimately be estahlished again. The time for 
this diffusive process to reach seme prescribed degree of ccanpleticoi 
will, of course, depend on the distance 8 through which it must act 
and on the coefficient of diffusion V of the fluid (kinematic 
viscosity for dynamical changes) . In particular, one-dimensional 
diffusion theory indicates that the time req.ulred is proportional to 
8 >. 

Temporal changes at the edge of a boundary layer will be evened 
out, or diffused, through the layer by this mechanism. Accordingly, 
the time S^/u will be taken as characteristic of the time required to 
diffuse throughout the layer a unit proportional change occurring at 
the outer edge. If the boundary layer is substantially the quasi-steady 
flow over a flat plate, 8^ is approximately proportional to '>x/u, and 
the "diffusion time" is proportional to l/u. 

Owing to the unsteadiness of the flow, changes in U appear at 
the edge of the boundary layer. The time required for a unit propor- 
tional change in stream velocity to take place cannot be given by a 
single quantity, but, arather, will be characterized by the following 
set of times: 



Accordingly, the ratio of "diffusion time" z/l7 to the time required 
for a unit proportional change to be Imposed is characterized by the 
following group of quantities; 

XU' aj' «y/^ 

These quantities are equivalent to the quantities (lO) , which may 
thus be considered a measure of the preraptness with which the boundary 
layer responds to impressed changes. Clearly, then, if the quan- 
tities (10) are very small, quasi-steady flow is to be e;xpected. If 
the velocity is increasing, the boundary layer becomes progressively 
thinner and hence responds more auad more quickly to changes at its 
outer edge. This effect clearly tends to establish quasi-steady flow. 
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Solution for Insulated Plate anfl 
Arbitrary Telooity-Tlme Eelatlon 

Solutions of tbe com.pressl'ble-flow equ.atlons (6) and (7) will be 
sought for the case of a flat Insulated plate traveling with a speed 
Tdiicb. my vary wltb time In a differentiable tbou^ otherwise arbitrary 
way. 

If consideration is restricted to a stage of motion where the 
diffusion time is already rather smll, the stream function may be 
defined as 




^ • • • ) 


T - *0^*1^»2^ . . . y vxx; 

1 Y Fu 

where O a — and is a function of X,T, the definition of 

which will be deduced in the course of the analysis. ' Similarly, 

1 r(0,So,Cp^;2; • • 0 (12) 

(The "recovery factor" is r(0,{;Q,Sp . . .),) 

Substitution of definitions (ll) and (12) into eq.uatlons (6) and 
(7) yields 


7TT* 




g» OP 

sa<>(j y Vbi + ^ i ^ ^ 


vq(j + Pr frj + 


i Pr(foo)2 . ^2 ^ r + i ^ Oro + 

Os 03 00 

U X/ 2 ^ X/ " 2 ^0 X/ 

n=0 n=0 n=0 * ' 


foC-^Cn) = 2; f<,(0,Cn) = f (O.Cn) » 0 
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Since the case of- an Insulated plate is to he treated, 
function r thus satisfies the boundary conditions 



the 




Eq.uatlons (13) and (14) may he omsidered correct only if the 
q.uantity V may he assumed constant. For constant wall temperature 
or for incompressihle flow, eq.uation (4) indicates that 0, and 
hence W, would he constant. However, under the condition of no heat 
transfer, the wall temperature, and hence W, must depend to some 
degree on both X and T. The circumstances under which this 
dependence may he neglected are discussed in appendix B. 

Equations (13) and (14) may he made self-consistent (that is, 
containing only functions of O and the various C n) 


Co® 




C2 = ^ > 


(17) 


These definitions correspond to the quantities (lO) previously derived 
by dimensional reasoning. If the quantities (17) are considered small. 


the functions f and r may be expanded as follows; 

f(0,5n) S '^(O) + + Clfl(o) + . * . + + . . . + 

CqCi^oI^^^ + • • • (^8) 

+ CifiC®) + . • . +Co^^Oo(^^ + . . . + 
CoCifoi^®) + • • • 

Substitution of equations (17), (18), and (19) into equations (13) 
through (16) and collecting terms multiplied by the various powers and 
products of the Jn yields, in part, 

pill + > =0 (20a) 

fo' ' * + FfQ' * - 2F*fo' + 3F» ‘fg = -4(2 - F*) + aOF* ' (20h) 

fjL* " + Ff^* ' - 4F‘fi* + 5F“fi = 4fo» (20c) 
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F‘(o.) B 2; 

F‘(0) = F(0) » 0 

(21a) 

fn‘H = ^n 

• (0) = fj^(O) = 0 

(21h) 

R«« + Pr FE' + ipr(F'*)2 « 0 


(22a) 

ro“ + Pr Ftq' - 2Pr F*rQ = Pr 

^-F*»fo»* + 8E + 2aE* - 3E*fo^ 

(22h) 

r]_ ' * + Pr Fr]_ ' - 4Pr F ‘r^ « Pr 

• 

^-F*»fi»» + 4ro - SE'f^^ 

(22o) 

• 

• 

E‘(0) = E(») = 0 


(23a) 

rn’(O) = rj^(o») = 0 


(23h) 


Eq.uatlcais (20a), (21a), (22a), and (23a) are tiie eo[uations 
appropriate to steady flow. Accordingly, eq.uatlons (18) and (19) indi- 
cate tliat for Htnn. 1 1 values of tiie the iamlnar compresel'ble boundary 

layer on an .insulated flat plate is nearly q.uasl -steady, ■with, respect 
to both Telocity and temperature profiles, provided that f^^ and r^^^ 

and their derivatives are of unit order of magnitude. 

Magnitudes of Cn* - The system, of eq,uatlons Just derived will he 
valid only under circumstances in which the are small, and (to he 

conservative) can he arranged in a decreasing seg.uence. The first, and 

XU* 

usually dominant member of this group, is to = — s-, Zh practical 

situations, the req.uirement of small ‘wlH commonly he met. For 

example, at a distance of 10 feet aft of the leading edge of a body 
having an Instantaneous acceleration of 100 feet per second per second 
(about three times the acceleration due to gravity), the q.uantlty {q 

will he less than 10 percent for a velocity greater than 100 feet per 
second and less than 1 percent for velocities greater than 316 feet per 
second. For the case of constant acceleration at 100 feet per second 
per second, £q will he less than 1 percent after the first 

3.16 seconds of flight. 

For constant acceleration, only nonzero £^. For 

fli^t at a velocity U proportional to []?“■, 
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and so fortli. For tMs case, the Cn taken in order clearly form, a 
decreasing seq.uenoe if Co itself is small and m is not too small. 
Since the successiye Cn ih-volve increasing powers of yelocity in the 
denominator, it will he true that the will generally form a 

diminishing seq.uence for moderate or high speeds, proYided that TI(T) 
is a different iahle function. The hl^er the speed, the more sharply 
U(T) may vary without violating this requirement. 

Solution of momentum equations. - The solution of equation (20a) 
subject to boundary conditims (21a) is available in reference 5. The 
following features of the function J(0) are Important: 

F“(0) « 1.328 (24) 


- 20 ] = - 1-721 


(25) 


Equations (20b) and (20c) subject to boundary conditions (21b) have been 
solved by direct numerical integration. The method used is described in 
appendix 0. Talues of foW. fiW. and their derivatives are presented 

in table I. The functions F‘(o), fQ*(a), and fi’(o) are plotted in 


figure 2. 

The wall shear stress 



following dimensionless form: 


may be obtained in the 


'w 





F"(0) +^fo”(0) +^2^fl“(0) + . . . 

U2 




or, according to equation (24) and the results presented in table I, 

T 


'w 




(0.664) : 


1 + (2.555) SI - (1.414) + . . . 


Tj2 


U'- 


(26) 


The leading teim of the right side is the quasi-steady solution. As 
would be expected, positive acceleration provides a skin friction higher 
than the quasi-steady value. 

Displacement thickness. - The displacement thickness S* may be 
defined as follows: 
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or, by use of eq.uatlons (s) and (ll). 



or, wlien equation (12) Is introduced. 



Wlien eq.uatlons (17) and (19) are used, tbe following eq.uation results: 



When eq.uation (25) and the information in tables I and II are applied 
for air bavlng a Prandtl number of 0.72, this eapresslon becomes 


1 + (0.645)(7-1) _ (q. 592) ^ [l+ (2.938)(7-l) mJJ- 

(0.694) [l - (2.702) (7-1) mJ] + . . .^ 


the leading term of which is the q,uasl-steady result. PosltiTe 
acceleration results in a boundary layer which is thinner than the 
q.uasl-steady boundary layer. 


Solution of energy eq.uatlone. - The solution of eq,uation (22a) is 
available ’(for esaniple, in reference 3), The remaining eq.uatlons (22b) 
and (22c), subject to boundary conditions (23b), have been solved for 
Pr = 0.72 by direct numerical Integration. The method used is described 
in appendix 0. As a check on numerical procedure, equation (22a) has 
been solved by the same method. Talues of P(0), rQ(<J), r2^(flt), and 
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tlieir first derivatives are presented in tatle II. The result 
R(0) = 0.848 agrees closely with the value 0.845 given in 
reference 3. The functions H(a), ro(o)^ and r]_(CJ) eire plotted in 

figure 3. The recovery factor for Pr = 0.72 is 

r(0) = E(0) + 23Z1 ro(0) + r^Co) + . . . 

or, from table II, 

r(0) = (0.848) 1 - (2.852) (3.266) 

I U'^ 

the leading term, of which is the q.uasi-steady solution 
Positive acceleration results in a surface temperature 
q.uasi-steady value. 


UNSTEADY PLOW IN PRESSURE GRALISrJT 

The foregoing analysis has shown that- the boundary layer on a flat 
plate traveling with a speed which varies with tlme_ becomes q.uasi- 

steady for vanishing values of- the parameters XU*/0^, X%’ */u^, and 
so forth. When theas_ q.uantitles are small but not zero, proportional 
correction terms arise. 

When the behavior of a boundary layer in unsteady flow with 
pressure gradient is to be investigated, it would appear expedient 
first to inq.ulre whether the boundary layer may be considered quasi- 
steady. If this is not done, needless effort may be expended on an 
essentially simple problem. In order to establish a criterion for 
quasi-steady flow in a pressure gradient, Jsnowledge is required of the 
significant parameters arising in the unsteady equations of motion. 


+ • 


» 


(28) 


(see reference 3). 
lower than the 


Integrated Momentum Equation 

The incompressible laminar boundary-layer equations for unsteady 


flow. 


Uf; + UU^J- + VUy = Uj. + TJUj. + VUyy 


(29) 


U^ + Vy = 0 


(30) 
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may te Integrated to yield 


^ (J^*) + ^ (u28**) + TJUj. 5* = V (31) 


where x and y are measured along and normal to the wall, 
respectively (fig. 1(c)), and where 

08 


6 *s 


08 


Jo 


(32) 


where 8 is considered a measure of the thlolmess of the boundary 
layer, as in the Pohlhausen analysis (see paragraph 60, reference l). 
The velocity profile is written 


^ s h(x,n)j 1 ) H 


(33) 


The differential eq.uatlon (29) may be evaluated at y = 0 to yield 

Ut +^x = -^(^)w 
whence, by use of definition (33), 

8% 8^+ 

-hr^^(x,0) = — :t- + 


V 


VTJ 


or 


where 


h^^(x,o) = -T 5 -(A + e) 


(34) 


8% 


A = ■■■ ("Pohlhausen A", see paragraph 60, 
reference l) 


I = 


8%J+ 


- VU 

lq.uation (34) indicates that the functional form, of h may be written 


§ = h(r,T]) 


(35) 
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whence , 


8* = 5H*(r); S** = 

Substituting eq.uatlons (35) and (36) Into eq.uatlon (31) yields 


/i * p dH*N a 

(2 ® arj sE-wj 


f 


tt 




e 




2 ) dS* 

dr 


th* + 1^ ^ = h^(r,o) 


(36) 


(37) 


Biuation (37) shows that, Irrespeotlye of the time variation of 
U “and 8, the boundary layer will be q,uasl-steady if the following’ 
quantities are each sufficiently small: 


S^t a 8%^ - 8%^.^. 

VU ^ ^2 " 


(38) 


As has been mentioned previously, the quantity 8^/v may be regarded 
as a measure of the time required to diffuse a change through the 
boundary layer. Changes in time of either U or 8 may give rise 
to effects to be diffused. Thus, the first and last of the 
quantities (38) may be regarded as the ratios of diffusion times to 
the times for unit changes of stream velocity to occur at a given 
value of 2 : because of acceleration and rate of increase of accelera- 
tion, respectively. The second member of the group is proportional to 


82 1 38 
8 St 


and thus relates diffusion time to time for the occurrence 


of a unit change In 8. 


Relative to a mass of fluid passing through the boundary layer, 
temporal-changes may occur by reason of variation of U or 8 with x. 
Measured in this way, the time for a unit change in velocity to occur 
is l/Uyr, I'ormlng the usual ratio with diffusion time yields 

s2 

which is the ’^ohlhausen A”. If this parameter vanishes, the 

boundary layer is that appropriate to the case of zero pressure gradient. 
It is interesting to note the analogy with unsteady motion: "Quasi- 
flat-plate " flow will occur even if U^c not vanish, provided that 

82 

Ux is small enough; the requirement is that changes in U diffuse 
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to the -wall umoiL more rapidly than they occur In the stream. Accordingly, 
the third member of the group (38) is analogous to the fourth hut 
depends on both Eulerlan and lagrangian time variations of Do 


It may further be noted that the quantity T (eoLuatlon (34)) is 
related to the diffusion throu^ the boundary layer of changes in shear 
stress due to stream velocity variation: Roughly, the shear is propor- 

tional to u/B, and because of velocity variation alone, its total time 
derivative is U{./S + UCJ2./S. The ratio of diffusion time 8^/u to the 

time for unit proportional change of shear to occur is thus 


q 2 m^\ 

» U\ 8 ~T~ ) " 


8 %^ 8%^ 

PU V “ ^ 




Ror the case of zero pressure gradient, the parameters and 

( "irj reduce to the previously discussed g.uantlty — g-, and the 

parameter — 5 — reduces to Therefore, it is to be eapeoted 


that 


8 %. 




PU 


t a 




■will usually be of the same order and, for 


moderate or hl^ speeds, ■will be much larger than either 5— or 

P^ 

p2u 


The foregoing discussion pertains directly only to incompressible 
flov althou^ similar conslderatlohs mi^t be expected to apply in 
compressible flow. 

Ror small, but nonzero, values of and 

proportional correction terms of the same order of magnitude arise in 
formulas (26) to. (28) for shear stress, displacement thictoiess, and 
recovery fac^tor in flow ■with zero pressure gradient. Since xU'/Cr 
and x%' */u^ special cases of the q.uantltles (38), it is 

expected that for cases involving pressure gradient, the magnitudes 
of the q.uantities (38) Indicate q.uantitatively the order of magnitude 
of departure from the q.uasl-steady condition. 

In a specific problem, the q.uantities (38) may be estimated as 
follows: Since a criterion for q,uasl-steadlness is desired, the 

boundary- layer thickness 8 may be estimated on the basis of q.uasl- 
steady flow. If the resulting -mlues of g.uantitles (38) are of such a 
magnitude as to Indicate nearly q.uasi-steady flow, the assumption made 
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in estimating 8 is proper. In many cases, this estimate of 8, and 
particularly of S8/3t, wMoli appears In the second memher of the 
group (38), can te made only hy finding the complete q.uasi-steady 
solution, a gamhle justified In vlerw of the great difficulty In finding 
a complete unsteady^ solution free of such assumptions. 


Oscillating Stream Telocity 


As an example to which the foregoing discussion applies, the 
stream velocity over a flat plate is represented as 


U = U]_ + Vli sin 2k k - tj 


(39) 


which denotes a flow wherein a sinusoidal velocity disturbance of 
amplitude and frequency k Is carried downstream with the mean 

velocity U]|_; the maximum magnitudes of the q.uantitles (38) are, from 
eq.uatlon (39), 


l8%+ 


VXS 


- UiB 


max 


2- 2 k Eeg ~ ~ ; Eb5 s 


B^xt 



p2 

max 



The quantity 


(^) 


max 


^2jt Eeg ~ j 


Un 


j cannot he estimated properly In advance, 

Tentatively, it may he assumed that it is of the_same order as 
8%t/^^ and then. If It is concluded that the magnitudes of the 
remaining g.uantlties indicate quasi-steadiness, the quasi-steady 
solution may he found and then analyzed to provide, as a check, a 

more precise estimate of 

Sir ( &k\^ ^1 

Thus, if Eeg ~ ^ « 1 and ^Eeg « 1, the flow in 

question may he regarded as quasi-steady. If, as a rough approximation, 
the formula for boundary- layer thickness for steady uniform laminar 
flow over a flat plate is used, and if it- is assumed that k = 10 cycles 
per second, TJ]_ * 100 feet per second, and uT^AJi = 0.1, then for 

approximately the first foot downstream of the leading edge , 

8k ^1 < 




UlUi 


0.01 


2229 



IfACA OT 2471 


19 


and 



< 0.001 

^1 


UHOTEADY TUEBDLEiriT B0DHl!Ai?r LAYER 


Tlie pliyslcal interpretation of tlie goyemlng paraaneters in is.Tnina.-i" 
flow would seem to have sufflolent generality to permit extension to the 
case of tlie turtulent "boundary layer. The most important quantities 
TOuld te 


. SX ■ §5 1 5 

V ^ uu ' 6 


(40) 


wiiere U is to "be Interpreted as an average diffusion ooefflolent. 

The coefficient of diffusion hy turtulent mixing may "be written Zu’, 
wiiere Z is a ^mixing lengbii” and is related to the Intensity of 

turbulent fluctuations. If the assumption is made that Z oc 8 and 
u* oc U, then the parameters (40) may he written 

^ ITx 8Ut &t 
^ U " jj2 ^ W 

^x 

In reference 6 it was shown that the parameter 8 may not he 

considered to govern the response of a turbulent boundary layer to an 
adverse pressure gradient. The evidence is incancluslve for favorable 
pressure gradient. Thus, the assumption made as to the form of the 
turbulent diffusion ooefflolent is of questionable value and only 
qualitative significance should be attached to the unsteady flow 
parameters 8 IT.j./D'^ and 8.(./o. In the absence of pressure gradient, 

8 may be taken to be (very roughly) proportional to xj thus 


8Ut 



00 




Therefore, it would seem that for zero pressure gradient, if 
is small enou^, a boundary layer will tend to be quasi- steady 
whether it is laminar or turbulent. 

Actually, for zero pressure greidlent, the quantity xIT^j/U^ may be 

obtained without assuming a specific form for the coefficient of 
turbulent diffusion. The boundary-layer thickness 8 may be thought 
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of as tlie result of lateral diffusion of Torticity from a point source 
(the leading edge of the plate) moving with a velocity TJ. (Compare 
paragraph 17 , reference 1.) If this process occurs without the com- 
plioating effects of pressure gradient, one-dimensional diffusion theory 

would indicate that 6^ = where V is either a laminar (molecular) 

or turbulent diffusion coefficient, and i/u is the time elapsed since the 

source of vortlclty passed the station x. Thus, the "diffusion time" 

b^/v appearing in the q.uantitles (39) would be given roughly by x/U, 

S^t -- uZ 

and the quantity reduces directly to • 

Of course, the parameter xTJ^/[J^ yields a rather unsatisfactory 

criterion for q.uasi-steadiness in the turbulent case because no indica- 
tion has been given of the magnitude of departure from q,uasi- steadiness 
associated with small nonzero values of this parameter. 


COHGLCJSIONS 

The laminar compressible boundary- layer flow over. a flat plate 
moving with a time-dependent velocity has been analyzed in detail. 

In the course of this analysis, a group of parameters arise, the 
magnitudes of which determine the nature of the flew unsteadiness. If 
the parameters are very large, the classical "start Ing-from-rest" solu- 
tion applies; if very small, the flow may be regarded as nearly q.uasl- 
steady; that is, at any j istant the motion is nearly that which would 
be obtained in steady flow at the conditions prevailing at that instant. 
The deviations of velocity and temperature profiles (for an insulated 
plate) from the q.uasl-steady state have been computed. Eelative to 
q.uasl-steady flow, constant acceleration results in a thinner boundary 
layer with greater skin friction and lower wall temperature. 

These parameters governing unsteady flat -plate flow may be shown 
to be special cases of a group of quantities which compare the time 
required to diffuse a unit proportional change in some stream quantity 
to the wall and the time for a unit proportional change in that quantity 
to appear in the outer flow. Eor moderate or hl^ speeds and rather 
large accelerations, the first member of this group is commonly very 
small. The generalized unsteady flow parameters appear eapllcitly in 
the Integrated momentum equation for boundary-layer flow with an 
unsteady pressure ^adient. 

In various problems of unsteady boundary- layer motion, for example, 
in the presence of oscillations in stream velocity or pressure, these 
parameters provide an advance criterion as to whether the flow need be 
considered essentially unsteady or whether the motion may be considered 
quasi-steady. 
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In. tlie atsence of a pressure gracLlent, the unstead 7 turhulent 
■boundary layer is espeoted to "be goTemed 'by the same unsteady flow 
parameter as gorems the corresponding laminar boundary layer. 

Lewis Fli^t Propulsicn laboratory, 

National Advisory Committee for Aeronautics, 

Cleveland, Ohio, May 25, 1951. 
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APFENDIX A - NOTATION 

The following notation is used in this report- 
A acceleration 

C constant appearing in temperature- viscosity relation 

c specific heat at constant pressure 

Jr 

F related to stream function for flat plate in steady flow 

f,g functions (identified by numerical subscript) related to 

stream function for unsteady flat plate- flow 

H*, H** dimensionless forms of S* and &** 

h dimensionless velocity profile 

M Mach nijmber 

Pr Prandtl number 

R function related to temperature profile for steady flat-plate 

flow 

r functions (identified by numerical subscript) related to temper- 

ature profile for unsteady flat-plate flow 

T, t time 

U stream velocity in x-direction 

u velocity in X, x-directlon 

V velocity in Y, y-dlrection 

X^ X coordinate along wall 

1) j coordinate nomnal to wall 

Z dimensionless coordinate (Y/^/vt) 

r* dimensionless parameter (A + 4) 

T ratio of specific heats 

6 boundary- layer thicloiess 
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5*, displacement and momentum thicknesses 

^ dimensionless coordinate ,(x/aT^) 

T) dimensionless coordinate (y/Q) 

9 temperature 


A 

|J. 

U 

e 

p 

a 




Pohlhausen parameter 
coefficient of viscosity 
constant ^ 




dimensionless parameter 
density 

dimensionless coordinate 
stream function 
Suh scripts; 

V evaluation at -wall (y = O) 

«*» evaluation in stream (y*^®») 




t 

vu J 




(■*^) 


Subscript notation for partial differentiation is used -when 
convenient. Primes are used to denote ordinary differentiation. 
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APPENDIX B 


ASSUMPTION OF CONSTANT V 


Equations (20) and (22) are derived for constant U . Deriva- 
tives of V with respect to time appear in the following terms of 
the differential equation (s) when definitions (ll) and (l2) are 
introduced; 





The quantity Vt/u * Crj/C may he obtained from equation (4). 

Since the flat plate is considered to be insulated and since this 
discussion is concerned only with orders of magnitude, may be 

taken equal to the stream stagnation temperature 

When it is noted that y for air is about 7/5, 



% Ct 1 u' 
^ = C = " 5 U 


Mod^ 

1 + i Mob^ 
5 


1 + i Mod^ - 


216° 


1 


216° R 
9 cd 


Thus, U|p may be neglected in expressions (bi) if 


_U *^T 
U' " 


i M<»2 


1 + 1 

5 6» 


1 + i MoF + 


« 1 


(B2) 


If = 0, 
is about 430° R, 


iL 

U* C 

jL 

U' C 


= 0^ and if Moe= », 


U ^tI 


1. If 00 


|U« C| 

= 0.07 for Moo» 1, and 0.25 when M,* “ 2. 


The errors in determination of fQ and Tq due to neglect of V 


are of order 


JJ % 
U* T 


T 

Therefore, for reasonable values of 6^ and 


for subsonic Mach numbers, the dependence of p on T may properly 
be neglected. Of course, for high speeds when ml gbt. be important. 
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ttie correction to the quasi-steady solution tends to he of slight 
importance since the magnitude of the correction is governed by such 
quantities as XU’/U^, vhich tends to be small for large velocities. 

Derivatives of y with respect to X appear in various terms 
of equations (s) -when definitions (ll) and (l2) are Introduced. 

The quantity appears in all cases in the factor 1 + X 

and may be neglected if X = X^ «1. Equation (4) combined with 

1? c 

equations (l2) and (l9) yields^ for small , 





1 S/" XD-»\ 
2TJ'c\^2 ] 


+ . 


If 1 * 0 ( 0 ) is taken to be of unit order, it may be shown that 
neglecting also involves an error in fQ and rQ of 

JL y 

U‘ C ’ 


order 
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APPENDIX C 
INTEGRAOnONS 

Integration of Mbmenttim- Equations 

The ordinary differential equation (20h) of third order for 
fQ(cr) and the houndaiy conditions (21h) constitute a two-point 

boundary- value problem, which may be solved numerically as follows; 

The scale of a is divided into intervals of 0.1. If the value of 
fQ(a) and its derivatives is known at five successive points, a 

fifth-degree polynomial may be matched to the five known and a sixth 
unknown value of f^" ' (a) . This polynomial may be integrated to 

yield fQ’’(o), fQ‘(a), and at the sixth point in terms of 

the six values of fQ‘''(a). The condition that fQ(cr) and its 

derivatives must satisfy the differential equation (20b) at the sixth 
point serves to determine at that point^the value of fQ*’'(tr), and 

hence the values of fQ’(cr), and fQ’*(a), Thus, given the 

solution for five consecutive points, the solution may be extended to 
the next point. 

The foregoing numerical procedure may be used for single-point 
boundary- value problems, provided f(d) is completely known for the 
first five Intervals. In this case, the first- five values are found 
by expanding fo(d) and p(a) in a Taylor's series about cr a o. 

The difficulty presented by the fact that boundary conditions are 
given at two points (cr«0,a>) is overcome by splitting the original 

The 


(Cl) 

be 


two-point problem into two single-point problems, as follows; 
solution fQ(cr) is written _ 

fo(a) 5 Cfo^^^ 

where C is a constant. The quantity (cf) is taken to 

the solution of the single-point boundary- value problem 


' ' ' + ' ' - 2F'fQ^^^ ' + 3F' » 0 

^ ( 1 )” 


(0) » li fQ(i) ' (0) = f^Ci) (0) = 0 


(C2) 
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and f, 




the solution of 


4 ^) 


t t t 


+ Ff^C^) 


- 2F 

( 2 )‘« 


+ 3F'‘fQC^5 » _ 4(2 - F') + 


2aF' 


( 0 ) 


^0^^^ ' Co) - f, 


0^^^ Co) - 0 


(C3) 


Definitions (Cl) to (C3) provide that 

fo(c^) = CfoC^J (a) + foC^) (a) will satisfy differential equation (20h) 
and will yield, at a «« 0, 


fo"Co}-c; fo‘(o) « fo(o) - 0 


If, after fQC^^ and fQC^^ have heen determined, C is taken eq.ual 
■to - \ the houndaiy condition fQ*(oo) »= 0 will "be 


O’ — i 03 


0 


Ccr) 


satisfied, and the solution for fQ(cr) may be written 


^‘oC^) 


lim fp^^^*Cg) 
cr-^o3 (1)' 

C«^) 


C^) + Co) 


The number of Intervals fOr which it is necessary to compute fQ^^^ 

sJiti fo^ ^ is primarily determined by the accuracy desired for the 
constant C. 

The solution of equation (2Pc) is also obtained by the foregoing 
method. 


Integration of Energy Equations 

Equations (22a) to (22c) are solved by the same method as is 
used for the momentum equations, except that a second-degree (rather 
than fifth-degree) polynomial is matched to three successive values of 
r' '(a). This simplification is considered warranted by the fact that 
the energy equations are of lower order than the momentum equations. 
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TABLE I - FUNCTIONS ASSOCIATED WITH VELOCITY FBOFILE 


CT 

^0 

f * 

^0 

f ’ • 
^0 

^1 

^l’ 

fl- 

0 

0 

0 

3.394120 

0 

0 

-1.878370 

.2 

.057747 

.529447 

1.953432 

-.036758 

-.359945 

-1.653824 

.4 

.194689 

,802508 

.829535 

-.138740 

- . 642866 

-1.147047 

.6 

.365840 

.881837 

.013897 

-.286337 

-.813374 

-.557151 

.8 

.538543 

.827366 

-.512697 

-.456399 

-.869151 

-.018114 

1.0 

. 691664 

. 693876 

-.783065 

-.627640 

-.828805 

.396729 

1.2 

.814136 

.528061 

-.845153 

- . 783550 

- . 721068 

. 653481 

1.4 

.903266 

.365638 

-.760117 

-.913840 

-.577909 

. 753609 

1.6 

.962232 

.229349 

-.594847 

-1.014405 

-.428269 

. 724592 

1.8 

.995747 

.128992 

-.409600 

-1.086252 

-.293752 

. 610518 

2.0 

1. 014474 

.063862 

-.247975 

-1.133778 

-.186499 

.461844 

2.2 

1.023130 

.026808 

-.130397 

-1.162855 

-.109292 

.313684 

2.4 

1.026420 

.008672 

-.057777 

-1.179289 

-.059187 

.193292 

2.6 

1.027286 

.001372 

-.019894 

-1.187871 

-.029603 

.108468 

2.8 

1.019088 

-.000727 

-.003696 

-1.192011 

-.013669 

.055563 

3.0 

1.018913 

-.000829 

.001437 

-1.193852 

-.005818 

.026086 

3.2 

1.018786 

-.000433 

.002081 

-1.194603 

-.002261 

.011299 

3.4 

1.018711 

-.000377 

.001358 

-1.194883 

-.000770 

.004460 

3.6 

1.018660 

-.000175 

.000668 

-1.195151 

-.000210 

.001640 

3.8 

1.018637 

-.000082 

.000276 

-1.195168 

0 

.000600 

4.0 

1.018622 

-.000043 

.000088 

-1.195160 

.000090 

.000250 

4.2 

1.018617 

-.000034 

.000014 

-1.195130 

.000120 

.000160 

4.4 

1.018612 

-.000031 

-.000012 

-1.195100 

.000140 

.000140 
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a?AELE II - FUNCTIONS ASSOCIATED WITH TEMEERATURE FROPILE 


a 

R 

R' 

^0 

1 

r * 

^0 

^1 

-l' 

0 

0.848000 

0 

-2.418820 

0 

2.769860 

0 

.2 

.835307 

-.126789 

-2.382630 

.374156 

2.678376 

-.846230 

.4 

. 797491 

- . 250398 

-2.267210 

.779209 

2.455221 

-1.335023 

.6 

. 735902 

-.362717 

-2.073428 

1.147468 

2.160495 

-1.573702 

.8 

: 653919 

-.452272 

-1.814377 

1.423885 

1.837448 

-1.634264 

1.0 

.557309 

-•.507224 

r-1. 512653 

1.569377 

1.513939 

-1.587060 

1.2 

.453863 

-.519942 

-1.196740 

1.567789 

1.207085 

-1.472302 

1.4 

.352151 

-.490678 

-.894813 

1.431645 

.928167 

-1.309995 

1.6 

.259817 

-.428239 

- . 630473 

1.200147 

.685411 

-1.113120 

1.8 

. 182067 

-.347426 

-.417568 

.926319 

.484160 

-.897970 

2.0 

.121126 

-.262489 

-.259398 

.660089 

.326110 

- . 684630 

2.2 

.076520 

-.185622 

-.150731 

.435205 

. 208880 

-.492490 

2.4 

.045900 

-.123264 

-.081592 

.265658 

.126810 

-.334450 

2.6 

.026133 

-.077049 

-.040868 

.149928 

.072500 

-.214940 

2.8 

.014108 

-.045402 

-.018722 

.077734 

.038410 

-.131330 

3.0 

.007206 

-.025243 

-.007696 

.056708 

.019510 

-.061160 

3.2 

.003468 

-.013248 

-.002726 

015507 

.010000 

-.034920 

3.4 

.001555 

-.006565 

-.000742 

.005627 

.004880 

-.017880 

3.6 

.000631 

-.003073 

-.000088 

.001555 

.002340 

-.008470 

3.8 

.000209 

-.001358 

.000053 

.000132 

.001200 

-.CX)3640 

4.0 

.000027 

-.000567 

.000032 

-.000242 

.000720 

-.001320 


r H.C . : 

i 

-.10974 

00 

^0 

- 2.99236 

OB 

1 r^^dcr 
0 

■ 3.22928 
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(t) Coordinates fixed in plate. 



(c) Curvilinear coordinates; y = 0 on wall. 
Figure 1. - Coordinate systems used in analysis. 



r*(a), fQ'ia), f^^a) 
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Figure 2. - Functions associated witli velocity profile. 
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Figure 3. - Functions associated with temperature profile. 
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